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We discuss in this paper the behaviour of minimal models of conformal theory per- 
turbed by the operator $13 at the boundary. Using the RSOS restriction of the sine-Gordon 
model, adapted to the boundary problem, a series of boundary flows between different set 
of conformally invariant boundary conditions are described. Generalizing the "staircase" 
phenomenon discovered by Al. Zamolodchikov, we find that an analytic continuation of 
the boundary sinh-Gordon model provides a flow interpolating not only between all min- 
imal models in the bulk, but also between their possible conformal boundary conditions. 
In the particular case where the bulk sinh-Gordon coupling is turned to zero, we obtain a 
boundary roaming trajectory in the c = 1 theory, that interpolates between all the possible 
spin s Kondo models. 
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1. Introduction 



The problem of minimal models perturbed by the $13 operator, both in the bulk 
and at the boundary, has been considered by various groups Nevertheless, various 

aspects are far from clear: for instance, the physical meaning of the solutions of the 
boundary Yang Baxter equation that have been found, or the relation with the boundary 
sine-Gordon model via quantum group truncation, if any. 

This is not only of academic interest. Minimal models arise naturally in quantum 
impurity problems, when one considers, for instance, situations involving several channels 
through a conformal embedding approach. When the bulk is massless and the boundary 
still perturbed by $13, one gets a boundary flow in the minimal model, ie an interpolation 
between two different types of boundary conditions (boundary fixed points). What sort of 
boundary conditions can be related in that way in the corresponding quantum impurity 
problem, and what is the boundary free energy along the flow, are questions of crucial 
importance (for recent progress in that direction using a different approach, see [Q.) 

Consider the most general problem of bulk and boundary perturbation of minimal 
models. By following the standard arguments ||, one finds that the theory with hamilto- 
nian 

H = H CF T + g rfx$i 3 $l3 + #b$is(0), (1.1) 
J — 00 

is integrable for any choice of the coupling constants g,gs- Here, Hq ft is the fixed point 

hamiltonian of a minimal model with central charge c = 1 , 6 . n , and $13 is the field 

with dimension /113 = 

±0 771+ 1 

In the case of a bulk theory, and for g > (so the bulk is massive) it is well understood 
how the problem can be considered as a truncation of the sine-Gordon model |5|] at coupling 
8^ = rJ+i ■ The simplest way to think of this is the scattering theory, where the truncation 
amounts to going from the soliton S-matrix (a vertex type solution of the Yang Baxter 
equation) to the kink S-matrix (an interacting round a face type solution) via quantum 
6j calculus based on the algebra sl q (2), q = exp (zvt 11 ^!). In this process, it is of crucial 
importance that the soliton S-matrix has a quantum group symmetry, in particular that 
it conserves the U(l) charge. 

It therefore appears very natural to try to describe the theory ( |1 . 1| ) with both a bulk 
and a boundary perturbation as a truncation of the boundary sine-Gordon model with 
hamiltonian 

H = H FB + G / da;cos^0 + GflCos^(0-0 o )(O). (1.2) 
J -00 2 
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Here Hps is the free boson hamiltonian. The coupling constants in (L2) are related with 
the ones in by G oc Jg. 

One quickly meets a major difficulty however: while the bulk S-matrix does conserve 
the i7(l) charge, the boundary R matrix, as worked out in ||, does not in generali. 
It is thus impossible to implement the quantum group truncation to obtain a boundary 
scattering for the kinks, and the question occurs, of whether the quantum group approach 
is any good to describe boundary perturbation of minimal models. 

The purpose of this note is to demonstrate that the problem (|1.1[) can be approached 
using sine-Gordon type models, but that to do so, one has to include boundary degrees of 
freedom of the Kondo type (see || for a discussion of the Kondo problem in the integrable 
field theory language). The idea is not so new, and was proposed independently and in a 
slightly different context by N. Warner in the case of supersymmetric models J7| . Some of 
the following results have appeared in very succint form in previous papers of some of the 
authors [H, U- 



2. Massless bulk 

We start by considering the case where the bulk is massless, ie g = G = 0. As is 
well known, this can be described using a massless scattering theory []10|] , []TI] . In the sine- 



Gordon case, one has massless solitons and antisolitons which are either R or L moving, 
with corresponding dispersion relation e = ±p = fie^ 3 , (3 the rapidity, \x an arbitrary energy 
scale, and scattering given by Srr = Sll = Ssg, the usual SG S -matrix, while the LR 
scattering is a simple constant phase. The truncated version of this is immediate: one 
has R or L moving kinks, with Srr = Sll = S m i n . Now, as far as the boundary goes, 
we claim first that action (|1 . 1| ) does not make sense without specifying the UV boundary 
conditions. We will restrict to conformal invariant boundary conditions (integrable flows 
only interpolate between such fixed points) . Then, a set of conformal boundary conditions 
was obtained in Jl2[ , |13|] . The simplest way to characterize them is to use the microscopic 



description of minimal models, where "height" variables take values on a A m Dynkin 
diagram, a = 1, . . . , m. The first type of boundary condition is obtained simply by fixing 



2 The only limiting cases where there is a conserved U(l) charge are the case of Dirichlet 
boundary conditions, Gb = oo, and the case of Neumann boundary conditions when the bulk is 
massless, G = Gb = 0. These correspond to boundary fixed points in minimal models, and do 
not describe the Bows one is looking for. 
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the heights to a constant; it corresponds, in the continuum limit, to a boundary state 
\h\ a > in the notations of |13|]. This state differs from \ hi jm+ i- a > by the phase of some of 
the coefficients only, as a result of the Z2 symmetry of the A m diagram. In the following, 
we will restrict to the case a < rn ^ 1 . The second type of boundary conditions is obtained 
by fixing the boundary heights to a value b and in addition, their neighbours one layer in 
from the boundary to a value c. For 6, c < setting d = inf(6, c), the corresponding 

boundary state is \hdi > (again, this state differs from \h m -d,i > by the phase of some 
coefficients only). For (6, c) = ^ + 1J and m even, the state is \h m / 2 .i > ■ In this 
paper, we will solve the problem of perturbing minimal models with the first type of 
boundary conditions. 

Our first claim is that the perturbation of the minimal model with UV boundary 
conditions of the type \hn > is, indeed, described by the only possible quantum group 
reduction of the boundary sine-Gordon model, that is, is trivial. A simple illustration of 
this is provided by the case m = 3, the Ising model. In that case, our boundary condition 
corresponds to fixed spins, while the boundary operator $13 corresponds to the boundary 
magnetic field: clearly, perturbing fixed spins by a boundary field is a trivial operation. 

Consider now the UV boundary condition \hi a > (a < 22 ^-), and set a = 1 + j. We 
conjecture that the problem is described by the quantum group truncation of the following 
hamiltonian 

H = H FB + G B (s+ e -*/W(0)/2 + g^i/W^/a j ? (21) 

where the opertors are in the spin j/2 representation of the quantum group sl q >(2), 
q' = exp (in^-^j , and = ;^pj- For Gb real, ie gs > 0, this problem is nothing but the 
anisotropic spin j/2 Kondo problem, which is well known to be integrable \TA\ [IS], [IB]. The 
R matrix is given by the Yang Baxter solution based on sl q {2) for scattering a spin 1/2 at 
rapidity (3 through a spin (j — 1)/2 at some fixed rapidity (3b (multiplied by the appropriate 
factors to ensure crossing and unitarity); the value of (3b is related with the coupling 
constant, [ie^ B oc g B 2 . Note the renormalization of the quantum group parameter and 
of the spin. The latter occurs because the R matrix description is really an IR one, and 
in the IR, the spin j/2 has been partially screened by the bulk degrees of freedom, so 
only the remainder (j — l)/2 appears. The case j/2 = 1/2 deserves special mention. 
In that case, there is no left over spin, and the R matrix is a simple CDD factor, R = 
-ztanh(£^ - f V 

The quantum group truncation of this R matrix is straightforward: it becomes a 
solution of Yang Baxter for scattering a kink with spin 1/2 adjacency rules through a kink 
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with spin (j — l)/2 adjacency rules. In the case j/2 = 1/2, we still get a CDD factor. To 
summarize 

R =Id, a = l 

# = -ztanhl — — l,a = 2 (2.2) 

R = J R 1 /2,(a-2)/2 (/3 _ 3 < a < m + 1 

The boundary free energy can easily be read off from the thermodynamic Bethe ansatz 
T7J]. As in the bulk case, we simply truncate the TBA for the bosonic theory - here the 



anisotropic Kondo model. Introduce the A m _2 Dynkin diagram with nodes labelled by 
an integer, and incidence matrix N pq ie N pq = 1 if the nodes p and q are connected, zero 
otherwise. Introduce the pseudo energies e p solutions of the system 

^^V-T^/f _ R l_ M ,n(l + e -^'>). ( 2 .3) 

Then one has, for the perturbation of the conformal boundary condition \h\ a >, with 
2 < a < and g B > 0, 



d/3 In (l + e- 6 "- 1 ^)) 
2^ cosh(/3 - p B ) 



til = -r / g „ (2.4) 



Of particular interest are the boundary entropies [0 in the UV and the IR, related to the 



degeneracy factors by s = Ing. From (|2.4| ), one has 

guv _ 1 + e -^-i(-°°) 
Slr ~ l + e-^-i(~) 

In the limit j3 — > ±oo, the e's go to constants. Setting x = e~ e ^~°°\ we have to solve the 
system 

x p = (l+a: p _ 1 ) 1 / 2 (l + x p+1 ) 1 / 2 

where xq = and x m -i = 0. The solution is well known 

' sin liE±!! x 2 
x p = I , m + 1 I - I . 



sm 



ra+1 



Similarly, for (3 = oo, we have to solve a similar system but with one less node, whose 
solution is thus 



7rp \ 2 

sin — \ 
1 sm — / 

m ' 



It follows that 



guv 
gm 



sm 



771+1 



sm — 

m 



sm m+1 sin -A— -2 



(2.5) 



On the other hand, recall the formula for the degeneracy factors of the boundary states 



9la 



.1/4 



m(m + 1) 



(-1) 



a+l 



sm — 
2 ^_ 



1/2 



7ra 



sin 



sin 



771+1 



m + 1 



and 



9di 



m(m + 1) 
We thus see that guv/gm = gia/g a - 



1/4 



("I) 



d+1 



sin 



777+1 



sm 



1/2 



7T(i 



sm 



rn 



i,i- 



From the microscopic interpretation given above, 



this is a flow from fixed boundary conditions on one row (a) to fixed boundary conditions 
(a, a — 1) on two neighbouring rows. Alternatively, by considering the row one layer in, we 
can think of it as a flow from a boundary condition where two heights a + 1 are allowed, 
to one where heights are fixed to the value a — 1. It is easy to check that the dimension 
of the UV perturbing operator is indeed his, since the diagram has the same skeleton as 
for the bulk perturbation. One also checks that the IR fixed point is approached along the 
direction <& 31 . 

When a > of course, one can still use the previous results by exploiting the 

symmetry of the lattice problem under a —>■ m + 1 — a. In microscopic terms, the flow is 
now from one height a fixed to two neighbouring heights fixed: a, a + 1. 

Nevertheless, the TBA diagram is not Z2 symmetric because the source term is only on 
the left. Therefore, if we use the same R matrices as before with a > m ^" 1 , we will describe 
another flow. By comparing with the boundary entropies, it is a flow from \hi a > (= 
|/ii jm _i_i_ a >) to I /id— 1,1 > (= |^m+i-a,i >)• In microscopic terms, setting a' = m + 1 — a, 
a' < m ^" 1 , we have a flow from one height fixed to a' to heights fixed a', a' + 1 

It is reasonable to expect that this flow corresponds to the perturbation with gs < 0, 
and this can be proven easily by considering the functional relations approach to the 
problem |IB[],||15||. Setting f^ r = — TlnZj, a = 1 + j, the Zj can be shown to obey fusion 
relations of the type 

Zi(GB)Zj(GB] 



B 



+ Z j+1 (q'-^G 



B, 



(2.6) 



For minimal models, these relations close, and one has 



z j(G 



Bj 



Z m -\-j\iGB) 



(2.7) 
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If we expand the boundary free energy in powers of the coupling constant gs oc G 2 B for 
either case, we get indeed identical expansions but with odd terms having the sign switched. 



To summarize, we have, for a < m ^" 1 , 



fbdA9B > 0, g = = - T / K — — — — , flow a^(a,a-l) 

J 2tt cosh(/3 - Pb) (2 g x 

i \ r dB In (l + e - e m- a (l3)\ 

f2(9B<0,g = 0)=-T / f- [ + J K flow«-(«,«+l) 

J 2tt cosh(/3 - p B ) 

where the e's solve the equations (|2.3|) . As in the bulk case, we expect the correspondence 
between the parameter Pb and gB to be the same in both regimes up to a sign, (ie^ B oc 

m + l 

\9b\ 2 • 



3. Massless flow in the bulk 

It is interesting to discuss the related problem where one has a massless flow both in 
the bulk and at the boundary. In that case, the bulk scattering is well known, and involves 
diffusion of left and right particles in a non trivial way We conjecture that the TBA 
is given by the following. Introduce the pseudo energies e p solutions of the system 

e P {P) = ile% 1 + ^5 p , m . 2 -Tj2N pq j g _ ^ In (l + e^) . (3.1) 



Then one has, for the perturbation of the conformal boundary condition \h\ a ) with 1 < 



a < 



f^l(9B>0,9<0) = ~ J 



dp 
2^ 



In (1 + e-^- 1 ^) In (1 + e" 6 — ^) 



cosh(/3 - p B ] 



cosh(P + p B ] 



(3.2) 



The mass scale \x now has a direct physical meaning. The other scale in the problem is 
the boundary scale Tb = j2e@ B . 

Various checks are available for this conjecture. First, consider the case where p, — > 0, 
so there actually is no flow in the bulk, or, more precisely, the whole boundary flow takes 
place at energies in the deep UV of the bulk flow. Then, we have to send Pb — > oo to 
keep Tb finite. In doing so, the TBA factorizes in two parts, corresponding respectively 
to the limits P — > ±oo. In the first case, setting Tb = ^e° B and fie 13 = \ie B , the mass term 
on the m — 2 th node of the TBA disappears, and we get for the boundary free energy a 
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contribution exactly equal to half the one in ( j2.4|) , with (3 b = Ob- The left moving part of 
the TBA gives the same result, so the two add up to what we had previously ( (pT4] )), as 
desired. 

Another interesting case is what happens at fixed (3b when jl is varied. As jl — > 0, the 
boundary is in the UV limit for the A m model, while as ft — > oo, it should be in the IR 
of the A m _i model. From our formulas, when jl — > 0, the two source terms in the TBA 
disappear, so we get, using results of the previous section 



fi 



bdr 



-Tin 



sin 



m+l 



sin 



(3.3) 



m+l 



On the other hand, as ft 
so 



oo, the two nodes of the TBA carrying source terms disappear, 

Tr(q-l) 



r (a) 
J bdr 



— T In 



sin 



m — l 



sin 



(3.4) 



TO— 1 



The difference of ( |3T3| ) and ( pT4] ) corresponds to a flow from \hi >a ) in the A m model to 
\h a -i^i) in the A m -i model, as expected. 

Finally, we can also consider the case where jl is very large and (3b — > — oo, so Tg 
is finite. Setting Tb = fJ-e 03 , by varying Tg we should describe the boundary flow in the 
A m _i model. Indeed, in that limit, for the first integral of ( [3.2| ), what contributes is the 
limit (3 — > — oo. Setting fi,e@ = fie d then, we get the same expression as (|2.4|) , while for the 
TBA the m — 2 th node has disappeared. A similar thing occurs for the second integral, so 
we get indeed the same thing as ( |2.4|) , but with one less node for the TBA, as expected. 



Of course, one can also check the limit where one is always in the UV of the boundary 
flow, or where one is always in the IR of this flow. 

The conjecture has an immediate extension to the case when grg < 0: 



fi2(9B<0,9<0) = ~ J f 



In (1 + e -e«-i(/9)) In (1 + e -^-a(/3)) 



cosh(/5 + /3 B ) 



cosh(/3 — (3b) 



(3.5) 



We leave it to the reader, to check that the various limiting cases are correctly reproduced. 



4. Boundary roaming trajectories. 

In |?D[ , Al. Zamolodchikov introduced a scattering theory with a single bosonic parti- 
cle (and no bound state) which appears to interpolate ( "roam" ) between successive minimal 
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models, exhibiting plateaux where the running central charge takes values c = 1 — ni ( m + , , 
as the RG scale is varied. 

This scattering theory can be formally considered as an analytic continuation to com- 
plex dimensions of the sinh- Gordon model 



H = H FB + u dxcosh(3(j). (4.1) 

J — oo 

Defining the coupling constant 

b 1 ? 



2ni + £i 
1T 



the continuation B — > 1 ± — 9n in the sinh-Gordon S matrix coincides indeed with the 



scattering theory in |2(J • Under this analytic continuation, each choice of sign leads to a 
different complex dimension A± for the perturbation in Q4.1|) : it was argued in [20] that 
the action corresponding to the continued scattering theory involves a real combination of 
both choices. 

Though the results of |20| are still to a large extent mysterious, it is interesting to see 
whether they extend to the boundary case. We thus consider the sinh-Gordon model ( f4.1| ) 
in the half space only, and add a boundary term 

# B = Acosh^(0(O)-0 o ). (4.2) 
The boundary reflection matrix is known for real (3 |2I] ; defining 

^ ' sinh(§ - 1 ' 

the reflection matrix is given by 

g(l)g(2-B/2)g(l + B/2) 
ai - E)i{\ + E)0l - F)i(l + F) y ■ ' 

with E and F related to the boundary parameters A and cpo in a complicated way. 

We now wish to formally extend the theory to complex values of f3, by continuing 
the coupling B in the bulk and boundary scattering matrices; for the reflection matrix to 
remain a pure phase, one needs to adjust the scales E, F appropriately in the complex 
plane. 

To compute the full boundary free energy, all terms of order 1/L (L the size of the 
system) have to be carefully taken into account. We are mostly interested in changes 
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of boundary entropy on the bulk plateaux, and for this, simpler formulas are sufficient. 
Introducing the bulk TBA 

-mR cosh 6 + e + ^-<f> * log(l + e" e ) = 0, (4.5) 

where the kernel is 

= -id e logS = — I H tt^ — ft, (4-6) 

cosh^y + t/ ) cosh(6/ — vq) 

one has 

Ss (E,F) = i- J ^ K (^)log(l + e- e ^), (4.7) 

where k = — ide log i?. The boundary flow at a fixed value of mR is completely described by 
the parameters E, F . For the rest of this work, we will further restrict to the case E = 0. 
corresponding to (fio = in the sinh-Gordon action, and the previous <pi3 considerations. 



In that case one has 



2 cosh cosh# 

K(9) = 5 2 - 5— =r- 4. 

V 7 (cosh 2 9 - sinh 2 2^ V 



To fix the ideas, let us start by reproducing one of the results of |2fJ: the function 
c(mR — r) is plotted for #o = 50 in Fig. 1. As can be seen, this function interpolates 
between the values for successive minimal models, starting from the Ising model with c = -| 
(m = 3). 

Let us now position ourselves on one of these plateaus by choosing r appropriately. 
We can then study what happens to the boundary entropy as we change F, which describes 
the scale of the boundary coupling constant A. The results are shown in Fig. 2 for the 
first three plateaus. 

A detailed study of the numerical values of the boundary entropies reveals that the 
boundary trajectory interpolates between all allowed flows in the minimal models discussed 
in section 2, though (roughly) half of these flows are "inverted" . Let us discuss this further. 
Consider first the case c = 1/2: the flow observed in Fig. 2 has a boundary entropy of 
log(2)/2 in the UV, and zero in the IR, corresponding to the flow from (in microscopic 
variables) 2 to 2, 1 in the As model (free to fixed boundary conditions in the usual Ising 
model). The next plateau shows more structure: the boundary entropy starts with a 
plateau at log(2)/2 and then rises to another plateau with value log (2 cos -|): the entropy 
difference is the opposite of the one for the flow from 2 to 2,3 in the A<± model: in the 
physical flow the boundary entropy would decrease, while here it increases. The running 
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entropy then decreases to zero, and the new entropy difference is the one for the flow from 
2 to 2, 1 in the A4 model. The next curve for c = 4/5 is getting even more intricate but 
each difference between the plateaus corresponds to one of the flows described in section 

2. From left to right, the difference of successive entropies correspond to flows from 3 to 

3, 2, the opposite of the flow from 2 to 2, 3 and finally the flow from 2 to 2, 1 in the A§ 
model. The same pattern have been verified for the next two plateaus of c, corresponding 
to m = 6, 7. 

These observations are somewhat satisfactory, bringing additional evidence that the 
roaming trajectory truly interpolates between minimal models. It is quite mysterious that 
the entropy actually increases along some of the flows, though this is not forbidden by the 
"g-theorem"p2|, since the theory we are dealing with involves complex dimension. The 
pattern is that, every time the sinh-Gordon model tries to "mimic" a flow of the type 
a — > a, a — 1 (a < in the A m model, it does it right, but every time it tries to mimic 

the flow a — > a, a + 1, it does it in reverse, producing a flow a, a + 1 — > a instead. 

It is especially interesting to consider the case m — > 00, ie when the bulk coupling 
vanishes. In that case, one expects a roaming within the c = 1 possible boundary condi- 
tions. Though this is difficult to study numerically, it is easy to find out what happens by 
taking the m — > 00 limit. For fixed a, as m — > 00, the ratio of g factors for the inverted 
flow a, a + 1 —>■ a goes to unity, while the ratio for the normal flow a —>■ a, a — 1 goes to 
The boundary entropy therefore interpolates between the logarithms of successive 
integers. Integer boundary degeneracies are well known in the c = 1 theory: they corre- 
spond to Kondo models with spin j, and g = 2j + 1. It is therefore tempting to conjecture 
that we have here a trajectory interpolating between successive spin j Kondo models, all 
the way from j = 00 in the deep UV (F — ► 0) to j = (F — > 00) in the IR 

Acknowledgments: we thank W. M. Koo and A. Leclair for useful discussions. 
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Figures Captions. 

Fig. 1: Roaming trajectory for 6> = 50. 

Fig. 2: Boundary flows for the three first plateaus of c(mR). 
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